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Copious production of top - anti top quark pairs at hadron colliders has enabled various probes into 
the properties and interactions of top quarks. Among the various presently measured observables, 
the forward-backward asymmetry (FBA) in ti production measured at the Tevatron significantly de- 
viates from the standard model predictions, and many models of new physics have been invented to 
explain the puzzle. We consider the consistency of the simplified single-resonance models containing 
a color octet axial- vector ( "axigluon" ) , color triplet or sextet weak singlet scalars, weak isodoublet 
scalar, flavor-changing neutral Z' , or charged W' vector boson with existing ti production mea- 
surements. Among the considered models only an axigluon can reproduce all Tevatron observables, 
without being in severe tension with the recent LHC results on ti production cross section, charge 
asymmetry and top-spin correlations. The LHC charge asymmetry measurements exclude the W 
and Z' explanations of the Tevatron FBA anomaly. On the other hand, all scalar models predict 
notable deviations in several top spin observables, and the recent top spin correlation measurement 
using the "helicity" spin quantization axis by ATLAS already provides a significant constraint on 
possible explanations of the Tevatron FBA anomaly. Future precise measurements of top spin cor- 
relations and especially top polarization could differentiate between scalar t-channel models, while 
they are less sensitive to pure axigluon contributions. 



I. INTRODUCTION 



Excellent performances of the Tevatron and recently the LHC have enriched our knowledge of top quark physics 
with many interesting results. Most intriguing are the observed deviations from the Standard Model (SM) predictions 
for the Forward-Backward Asymmetry (FBA) in ti production at Tevatron. At the same time the measured inclusive 
ti production cross section and its ti invariant mass (m t i) distribution are in a good agreement with SM predictions. 
Similarly, early ti production Charge Asymmetry (CA) measurements at the LHC show no deviations from zero, 
consistent with tiny SM predictions. This observed pattern of deviations motivated many studies of New Physics 
(NP) aiming to explain the puzzling phenomena. The required size of possible NP contributions to the FBA points 
towards tree-level effects and most studies employ the description in terms of s, t or it-channel resonance exchanges 
(c.f. pQ for a recent review). There are several candidate states which can serve to explain the Tevatron experimental 
results, such as the s-channel axial- vector color octet ("axigluon"), or alternatively u, or i-channel exchanged scalars 
(isosinglet color triplet or sextet, and color-neutral weak isodoublet) or flavor-changing neutral (Z') or charged (W) 
vector bosons. 

In order to shed light on the possible role of NP in top - anti top pair production at the Tevatron and the LHC, it 
is important to investigate all possible observables. For more than a decade proposals have been made to investigate 
top quark spin polarization and correlation effects. Namely, due to its very short lifetime, top quark spin information 
is not spoilt by hadronization effects and can be reconstructed from the angular distribution of decaying particles [2] . 
With the large ti datasets at the LHC, the possible presence of NP in ti production might be tested through new 
variables such as the (anti)top spin polarization and top-antitop spin-spin correlations [3H5]- The intricate relations 
between the ti charge asymmetric production, the underlying helicity amplitude structure and angular distributions 
of top quark decay products in ti production have also been pointed out recently in [6l [7] . 

Motivated by these developments and recent first precision measurements of CA and top spin observables at the 
LHC we reinvestigate effective single resonance NP models which aim to address the Tevatron anomalous FBA. 
Performing a fit of NP model parameters taking into account the presently measured ti production observables at the 
Tevatron and the LHC we determine the viable parameter regions in each model (previously, a similar analysis using 
the first LHC data has been performed in [§]). In this we pay special attention to the apparent tension between the 
large FBA values as measured at the Tevatron and the existing tight constraints on the CA from the LHC. For the 
viable model parameter regions we predict top spin observables at the Tevatron and the LHC, evaluate the impact of 
existing measurements, and point out which near future experimental analyses could discriminate between the various 
NP proposals for the FBA puzzle. 

This paper is organized as follows: in Sees. [H] and |III| we introduce the effective NP resonance models addressing 
the FBA anomaly, and the relevant Tevatron and LHC observables constraining these models, respectively. The 
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concordance of the NP models with existing tt production measurements is quantified in Sec. IV and the corresponding 
best fit regions are identified. In Sec. [v]we discuss the top spin phenomenology in hadronic tt production and give 
the FBA correlated NP model predictions for the relevant top spin observables at the Tevatron and the LHC. Finally 
we conclude in Sec. [VI] 



II. MODELS 



We consider NP models which affect the tt production at hadronic colliders at the tree level and consider a single NP 
amplitude interfering with the SM contributions at a time. Such scenarios can then be classified according to the new 
resonances coupling to quarks and exchanged in s— , t— or u— channel. We note however that our analysis is applicable 
also to NP models with more new degrees of freedom accessible at low energies (as for example in approximately 
flavor symmetric models [9]) but where the tt production phenomenology is dominated by the exchanges of a single 
intermediate state. Among the plethora of possible spin, weak isospin, charge and color assignments, only a few of 
such states can produce a sizable positive FBA at the Tevatron without being in gross conflict with the measurements 
of the total cross-section and/or the m t i spectrum [TO]. These include an s— channel exchanged (axial-) vector color 
octet (axigluon G') boson [TTHL3"] . neutral (Z') [T3] or charged (W) [15] vector bosons, coupling chirally to quarks 
and exchanged in the t— channel, a scalar isodoublet |I6j whose neutral component (</> ) contributes in the t— channel, 
as well as scalar color triplet (A) jTYJ, [TH] or sextet (E) [T7] coupling chirally to up-type quarks and contributing in 
the u— channel. The relevant interaction Lagrangians for all the considered models can be found in Appendix A. 



III. OBSERVABLES & METHODOLOGY 



We probe the parameter space of considered NP models with presently measured observables in tt production both 
at the Tevatron and the LHC. First, there is the inclusive FBA as measured by both the CDF [19] and D0 [20] collab- 
orations. Combining all the measurements on independent datasets by both experiments and adding the uncertainties 
in quadrature leads to a naive average 

A FB = 0.187 ±0.037, (1) 

compared to the NLO QCD prediction 19, 2TJ [22] (consistent with recent NLO+NNLL estimates [23] within errors) 
including leading electroweak (EW) contributions [24[ [25] of Apg — 0.07(2). At the LHC, the most precise mea- 
surements of the Charge Asymmetry (CA) to date by ATLAS [33 and CMS [57] can again be naively combined to 
yield 

A c = 0.001 ±0.014, (2) 

and compared to the SM prediction of A C M = 0.007(1) [5TJ [53 [MIS] ■ 

On the other hand, while both CDF and D0 have also presented FBA measurements in bins of m t j and top - anti 
top rapidity differences, only CDF [T9] presents results unfolded to the partonic (truth) level with the result for two 
m t i bins 

A' FB = A FB {m ti < 450 GeV) = 0.078 ± 0.054, (3a) 
A h FB = Apsirtiti > 450 GeV) = 0.296 ± 0.067, (3b) 

to be compared with SM predictions [H [52 [53 [5J [55J of Af™ = 0.05(1) and A h p B M = 0.11(2). Since these 
measurements are not independent of the inclusive FBA average, we take two complementary approaches. In the 
first (approach A) we use the inclusive A FB , while in the second (B) we instead consider only the two binned FBA 
values A FB (we take the associated uncertainties as uncorrelated). In both approaches we also treat the production 
cross-section measurements a bit differently (see below). We compare the constraints on NP model parameters in 
both approaches (A and B) in the next section. Similar m t i binned CA measurements have been performed by 
ATLAS [28] and CMS [27], but with large uncertainties, which make these observables at present less constraining 
than the inclusive CA measurements. 

Also crucially important observables at the Tevatron are the total inclusive tt production cross-section, with a 
recent experimental average of [59"] 



cr TEV = (7.5 ± 0.48) pb , 



(4) 
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compared to the recent NNLO+NNLL QCD prediction of cr^ v = 7.07(26) pb [3D]; as well as the m t i spectrum 
(dcrxEV / 'dmtt) measured in several bins by the CDF collaboration [31] and in good agreement with NLO+NNLL 
QCD predictions [35]. 1 Since again the two observables are not independent, we consider one of them in each of 
the two complementary approaches (A and B) discussed above. In the first (approach A) we use otev but also the 
next-to-last bin (m t j £ [700, 800] GeV) of the CDF m t i spectrum measurement [3T] 

4 EV = (80 ± 37) fb , (5) 

to be compared with the SM expectation of p^ev = 80(8) fb [35]. While the two observables are not completely 
independent, the effect of the high m t i tail of the spectrum on the total cross-section measurement is negligible (i..e 
much smaller than the experimental error of otev)- In the second approach (B) we do not use the total inclusive 
cross-section cttev at all, but instead employ all the m t i bins of the measured CDF spectrum. In this we consider the 
statistical uncertainties as uncorrelatcd among the bins, but we use a 100% correlation approximation for both the 
systematic as well as SM theoretical uncertainties. 

At the LHC, the most precise inclusive ti production cross-section determinations are presently provided by AT- 
LAS [33J and CMS [33] when combining their dilepton, single-lepton and all-hadronic final state analyses. Again 
performing a naive average of both experiments, we obtain 

(tlhc = (172 ± 10) pb , (6) 

in good agreement with the approximate NNLO QCD prediction of uf^ c — (163^q) pb [35]. On the other hand the 
existing high m t i spectrum measurements 36J arc limited by moderate statistics and significant systematic uncertain- 
ties |37| . In particular, in the m t t > 1 TeV region they yield the constraint [36 

4f§ < 2 - 6 @95%C.L.. (7) 

°LHC 

Since the dominant sources of uncertainty should be significantly reduced with more data, the m t j spectrum mea- 
surements could yield important constraints on NP contributions to ti production in the near future. 



Finally, ti spin correlations (as defined in Sec. VA| have been observed both at the Tevatron [38, 39 and the 
LHC [301 IM] • At the Tevatron presently the most precise determination has been performed in the "beamline" basis 
by the D0 collaboration [3j|], obtaining 

C™= 0.66 ±0.23, (8) 

in agreement with NLO QCD prediction of SM = 0.78^q'q4 [31] • On the other hand, the most recent "helicity" 

basis analysis by CDF [38] yields 

C£f v = -0.60 ± 0.52 , (9) 

where we have combined the statistical and systematic uncertainties in quadrature. Within the large uncertainties this 
is in agreement with the SM expectation of C™g M ~ —0.35 [42]. At the LHC, the recent ATLAS measurement [41] 
of the "helicity" axis spin correlation coefficient 

C£J° = 0M±° o f s , (10) 

in agreement with the SM prediction of Ct^ ' M ~ 0.31 [35] already provides a significant constraint on possible 
explanations of the Tevatron FBA anomaly. A detailed discussion of the impact of present and future measurements 
of tt spin observables on the NP models under consideration is given in section |Y| 

Some of the models we consider can also be constrained by phenomenology not directly related to ti production like 
same-sign and single top production, dijets, electroweak precision observables and rare B processes (c.f. [I] for a recent 
review and also |43] ) . Since in the present analysis we want to focus on the discriminating power of ti observables at 
the Tevatron and the LHC, we will assume that other constraints can be evaded in suitable UV completions of the 
considered effective single-resonance models. 



1 It has been pointed out recently, that sizable EW corrections can reduce the SM da jdm t i predictions in the high m t j region 1251 . thus 
leaving more room for potential NP contributions. However, such effects are smaller than the current experimental uncertainties [31| 
and can at present be safely neglected. 



4 



We compute the SM+NP contributions to the selected observables at the partonic level employing the MSTW2008 
parton distribution functions |44j at fixed factorization and renormalization scale /j,f = fJ-B. = rn t = 172.5 GeV - the 
top quark mass used also in the above quoted measurements to which we compare our predictions. For all observables 
we compute all interfering SM (QCD) and NP contributions at leading order in a s (at the tree level). We then 
normalize our SM+NP tree-level estimates to the pure SM result, i.e. we define normalized cross-sections a 



a 



LO 



CT QCD 



and compare these to the experimental cross-sections normalized to the state-of-the-art SM predictions, 



(12) 



In this way we hope to capture a (universal) part of the higher order QCD corrections to the observables under 
study. 2 In the case of forward-backward, charge asymmetries and top spin observables which are defined normalized 
to the production cross section, we first subtract our pure LO QCD estimates to obtain an estimate of the NP (and 
interference) contributions, i.e. for observable O wc define 

A0 th = lo _ lo d _ ( 13 ) 

We compare these to the experimental values for the asymmetries and top spin observables from which we have 
subtracted the state-of-the-art SM predictions, corrected for the change in the normalized cross-section 

AO exp = O cxp - 4v7Csm ■ (14) 

tj th 

In practice the strong experimental bound on NP contributions to cttev in 

Eq. @ restricts a th ~ 1 for inclusive 

observables (especially at the LHC due to the dominance of the gluon-fusion subprocess, not affected by NP in all 
models under consideration) and the associated corrections to AO exp are always small. 

A second issue regards experimental acceptance corrections due to limited rapidity coverage of the CDF and D0 
detectors. In deconvolving the measured detector level observables to the truth (partonic) level, the existing analyses 
assume (anti)top rapidity distributions resembling the SM predictions. However, NP models with light mediators 
exchanged in the t— and u— channels can exhibit a forward scattering peak thus enhancing (or suppressing) the 
cross-section relative to the SM at high rapidities [UJ . To correct for these effects we employ the procedure suggested 
in [9] by computing and applying efficiency corrections to the NP cross-section distributions in bins of top - anti top 
rapidity difference and m t j . 

In addition to the above mentioned acceptance effects, t— or u— channel models can contribute to ii+jet produc- 
tion for resonance masses above m t via associated resonance-top production. This could in principle affect also the it 
measurements by increasing the reconstructed cross-sections |47j . The magnitude of such effects however is propor- 
tional to the branching fraction for the decay of the resonances into t+jet which could in principle be reduced in UV 
completions of the effective models by the presence of other decay channels. In addition their relevance also depends 
somewhat on the details of the experimental analyses. Consequently we do not include them in our procedure. 

Finally, t— or u~ channel resonances with masses below mt can contribute new decay channels of the top quark, 
thus enhancing its width. While such effects are constrained by the single top production measurements |48| . the 
exact sensitivity of existing measurements is difficult to asses. Thus we do not strictly impose any such constraint, 
but do highlight such effects if significant in certain part of NP model parameter space. On the other hand, all of 
the present tt production experimental analyses rely on the kinematic reconstruction of top quarks assuming the 
standard decay mode t — > bW. In particular they require tight reconstruction criteria for the W in both hadronic 
and leptonic final states (c.f. 05] )• We will assume that any new decay mode of the top quark would not pass the 
top reconstruction criteria thus effectively reducing the measured cross-section. We take this effect into account by 
multiplying the predicted cross-sections by the square of the ratio of the total SM top width (assuming V t b = 1) and 



the NP enhanced top width (both computed at leading order in QCD). In particular, we correct Eq. (11) to 

r LO 



* th - 4o- Br (* m ? • ( l5 ) 

°QCD 



2 Presently, leading QCD corrections have only been computed within the W' 1451 and axigluon I4GI NP models and for a limited set of 
observables. These known results suggest that higher order QCD effects in NP contributions are indeed similar in size to SM higher 
order QCD corrections. 
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FIG. 1: ti production constraints on the axigluon model in approach B, showing "low" (left plot) and "high" (right plot) mass 
regions: binned FBA at la in thick full green line, inclusive CA at la (2a) in thick dashed green line (thin dashed red line), 
CDF m t f spectrum at 2a in thin red dotted line. Parameter regions where the model can improve the SM x 2 by -A X 2 > 0,1,4 
are shaded in red, yellow and green respectively. 



IV. MODEL FITS 



In constraining the NP models from tt production phenomenology we combine the observables discussed in the 
previous section into a global \ 2 fit (Afbj Ac, otev, °tevi ^lhc and Ol HC m a PP roacn A, or Ap B , Ac, olhc, 0-lhc 
and (ifiTEV / dm t i in approach B as explained in the previous section) . While we do not assign a statistical significance 
to the particular \ 2 values, we consider NP model parameter regions as acceptable if they improve upon the \ 2 
obtained in the pure SM (xa sm = 1-9/d.o.f. in approach A and x% sm = 1-6/d.o.f. in approach B). For these regions, 
which represent scenarios addressing the FBA discrepancy and being in reasonable agreement with other existing ti 
production measurements, we predict the relevant top spin observables at the Tevatron and the LHC in Sec. |Vj 



A. Axigluon 



Among the possible interactions of a heavy color-octet (axial) vector G , only the purely axial couplings to uu 
(gX) and tt (g^) currents contribute to the FBA at the tree-level. 3 We thus consider their product g\g\ together 
with the axigluon mass (mg) as the free parameters of the model. Since the effects of a s— channel resonance on 
the m t t spectrum depend crucially on its width we follow the approach of [121 113] and assign an axigluon width of 
^a/ m G ~ 0.2, deferring the explanation of such a large width to a UV completion of the model. The singular features 
of an s— channel resonance on ti spectrum also make the approach B (using the m t i spectrum instead of the total 
cross-section and inclusive FBA) much preferred for the axigluon model. The results of the global \ 2 fit together 
with the impact of individual constraints is presented in Fig. [T] We single out two interesting parameter plane regions 
as the "low" and "high" mass region focused on the left and the right plots respectively (consistent with results of 
previous studies, c.f. [50]). We observe that even with the large unexplained axigluon width, axigluon masses which 
fall within the measured CDF m t i spectrum bins are significantly constrained and in particular cannot improve upon 
the SM agreement with current experimental results. On the other hand the finite axigluon width effects become 



Maintaining purely axial quark-axigluon couplings in an elcctrowcak gauge invariant way requires identical couplings to down quarks. 
However, since at the Tevatron and the LHC, the corresponding dd parton lumionsities are much smaller than uu, we neglect the small 
dd initial state axigluon contributions to our observables. 



6 




100 150 200 250 300 350 400 450 100 150 200 250 300 

m z [GeV] m w [GeV] 

FIG. 2: ti production constraints on the Z' (left plot) and W' (right plot) models in approach A: FBA at la (2<r)in thin full 
green line (thick full red line), CA at la (2a) in thin dashed green line (thick dashed red line), cttev at 2a in thin red dotted 
line, Ot EV at 2a in thin red dash-dotted line. In no parameter region can either of the two models improve upon the SM fit 
due to the tension between the FBA and CA measurements. 

insignificant in the "low" region around and below ti production threshold (for mc < 450 GeV, consistent with the 
findings of [13]) and in the "high" region above the kinematic reach of the Tevatron (for tuq > 700 GeV, as pointed 
out in [H]). At least in the later, CTlhc ES] already contributes a relevant constraint for niQ > 1 TeV and future 
precise LHC ti spectrum measurements could have a significant impact. We also note that for this model the tension 
between FBA and CA measurements is apparent but can be reduced to the la level. The tension is expected to be 
larger when symmetric couplings to d quarks (as formally required by SM gauge invariance for a pure axigluon) are 
included [ST] , 

B. Z' and W 

For the Z' and W models we consider only right-handed flavor changing u — t and d — t couplings (denoted by 
fa) respectively since these are sufficient to induce a positive FBA contribution. The resulting parameter fits for 
approach A are represented in Fig. [2j We observe immediately that while a consistent description of ti observables 
at the Tevatron can be obtained (in accordance with previous findings [52 ), these models cannot improve upon the 
SM results once the LHC measurements are taken into account. This can be traced to the inherent tension between 
the FBA and C A measurements (and thus independent of approach A or B to the cross-section measurements) . For 
Z'(W) these are tightly correlated (as previously pointed out in [ID]), as also illustrated in Fig.|3]where we have varied 
fn within the perturbative regime (\Jr\ £ [0,47r]) and the Z' (W) masses in the region m z <(w') S [100,500] GeV. In 
particular from these plots it is clear that the tension cannot be reduced below the 2a level. We also note that the 
tension increases slightly towards higher Z' (W) masses (when marginalized over the couplings). Consequently, we 
conclude that these models are excluded as viable FBA explanations and do not consider them any further. 
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FIG. 3: Correlation between the FBA (AA FB ) and CA (AA C ) contributions of Z' (left plot) and W' (right plot) models in 
shaded narrow gray bands. The thickness of the bands is given by the Z' (W) mass variation in the range [100, 500] GeV. The 
Tevatron average of FBA measurements minus the SM prediction at la (2a) is represented by the vertical blue band (dashed 
vertical lines). The ATLAS measurement of CA minus the SM prediction (the SM reference point is marked with "*") at la 
(2a) is represented by the horizontal green band (dashed horizontal lines). The tension between the two observables increases 
with increasing mediator mass (inner edge of the band corresponds to lowest mediator mass). 



C. Scalar isodoublet 



In this model the parameter scan is performed over the ubXl coupling (j/13) and massfm^) of the neutral isodoublet 
component <f)° . 4 The resulting parameter fits for approach A and B are presented in Fig. El We observe that in approach 
A (using cttev and c^ EV ) there is a small preferred parameter region at low 4>° masses (m^ < m t , consistent with 
findings of [H>]). Interestingly, is bounded from below by the LHC inclusive cross-section measurement. This is 
because the new t — > <fi°u decay channel tends to suppress the reconstructed cross-sections (see Eq. (15)) irrespective 
of the production mechanism inducing a tension between the Tevatron and LHC results. When considering instead 
the complete CDF m t i spectrum (approach B), the < m t region is no longer favored. This happens because a 
light scalar exchanged in the t— channel tends to harden the da/dm t f spectrum compared to SM QCD predictions, 
inducing a tension with the binned spectrum measurement by CDF (in this particular case, the discrepancy is most 
pronounced in the lowest m t f bin which was measured precisely by CDF |31j). Given our rough treatment of the 
correlated uncertainties entering dtrxEV /dm t t and also the neglect of subleading SM contributions [25] , we shall keep 
the low mass region, disfavored by the da/dm t i observables, in our analysis. What remains in addition is a less 
interesting intermediate to heavy <fr mass region (m^ > 200 GeV) where the FBA cannot be reproduced at the 
la level however the overall agreement with the present experimental data is improved with respect to the SM (in 
order to clearly distinguish it from the light (f>° region, we shall only consider parameter values where the overall fit 
improvement with respect to the SM is at least —Ax 2 > 4, shaded green in right plot of Fig. [4J. Finally in this model 
the tension between the inclusive FBA and CA measurements can be reduced below the 2a level, even though the 
68% CL regions for both observables cannot be reached simultaneously. 



4 We do not consider the potential dd — ¥ tt contributions mediated by the charged isodoublet component, which are suppressed by smaller 
parton luminosities. 
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FIG. 4: ti production constraints on the scalar isodoublet model in approach A (left plot) and B (right plot): FBA (inclusive 
in left plot and binned in right plot) at la (2a) in thick full green line (thin full red line), inclusive CA at la (2a) in thick 
dashed green line (thin dashed red line). In left plot otev at la (2a) in thick green dotted line (thin red dotted line), ctlhc a t 
la (2a) in thick green dash-dotted line (thin red dash-dotted line). In right plot CDF m t i spectrum at la (2a) in thick green 
dotted line (thin red dotted line). Parameter regions where the model can improve the SM x 2 by -Ax 2 > 0, 1, 4 are shaded in 
red, yellow and green respectively. Finally in left plot for m$ < mj, the contours of constant branching fraction Br(£ — > u<j>°) 
are displayed in thin gray dashed lines. 

D. Scalar color triplet and sextet 

We consider chiral couplings of A and £ to right-handed di-quark currents with a coupling constant 1713 . The 

results of the corresponding parameter fits are presented in Figs. [5] and [6] Both of these models experience a tension 
between the FBA and the ti spectrum measurements |47j . in particular at high m t j. While in approach A an overall 
improvement in the goodness of fit can still be obtained in some parameter regions of the two models (see left plots in 
Figs. [5] and [6]); taking into account the complete binned m t i spectrum measurement by CDF removes these preferred 
regions and the FBA discrepancy can no longer be addressed. Again, taking into account the possible caveats of the 
do I dm t i constraints, we shall keep the A and £ model parameter regions preferred in approach A in our analysis of 
top spin observables. 

V. TOP SPIN OBSERVABLES 
A. Definitions 

Since the spin information of a decaying top quark is not diluted by hadronization it is possible to investigate 
various top-spin polarization effects at hadron colliders. The starting point is the calculation of the production spin 
density matrices for the dominant partonic production channels of top quark pairs, 

q{Pi) + q(P2) -> t(ki,s t ) + t{k 2 ,sj) , 

g(pi) + g(p2) -> t(h, s t ) + t(k 2 ,s T ) , (16) 

via the relevant scattering amplitudes Mj where I = qq, gg. The differential cross sections for these partonic processes 
can be written as daj(st,Sj) — $ i\ 2 dT ' t i where is the relevant initial state flux normalization and dT t ^ is the 
ti phase space differential. The problem is thus reduced to evaluating the absolute squares of the relevant polarized 
scattering amplitudes 

\Mi\ 2 = ^Tr[ PJ (l + s t • a) ® (1 + %• a)] , (17) 
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FIG. 5: ti production constraints on the scalar color triplet model in approach A (left plot) and B (right plot): FBA (inclusive 
in left plot and binned in right plot) at la (2a) in thick full green line (thin full red line), inclusive CA at la (2a) in thick 
dashed green line (thin dashed red line). In left plot otev at la (2a) in thick green dotted line (thin red dotted line), cr TEV at 
2a in thin red dash-dotted line. In right plot CDF m t i spectrum at la (2a) in thick green dotted line (thin red dotted line). 
Parameter regions where the model can improve the SM x 2 by -A X 2 > 0, 1 are shaded in red and yellow respectively. 
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FIG. 6: ti production constraints on the scalar color sextet model in approach A (left plot) and B (right plot): FBA (inclusive 
in left plot and binned in right plot) at la (2a) in thick full green line (thin full red line), inclusive CA at 2a in thin dashed 
red line. In left plot otev at la in thick green dotted line, CTtev at lc in thick green dash-dotted line. In right plot CDF m t i 
spectrum at 2a in thin red dotted line. Parameter regions where the model can improve the SM x 2 by -* 2 > 0, 1 are shaded 
in red and yellow respectively. 

where pi is the corresponding partonic production spin density matrix describing the production of (on-shell) top 
quark pairs in a specific spin configuration, while s t (§j) is the unit polarization vector of the top (anti top) quark in 
its rest frame and <x = (<Ti, 02, cr 3 ) T is a vector of Pauli matrices. 

Conveniently, one starts with the most general decomposition of the spin density matrix p (suppressing the initial 
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state subscript I) for the production of a top quark pair 

p = Al ® 1 + B\ en ® 1 + B\ 1 ® cr,- + Cy cr l ® (j j , (18) 

where the functions A, £?* (i?*) and C»j describe the spin-averaged production cross section, polarization of a top 
(anti top) quarks and the top - anti top spin-spin correlations, respectively. Using the spin four- vectors defined as 

k x • s t „ ki(ki • St) 



m t ' m t (E t +m t ) 

4 = (^A+ ^lgj , 



m t ' m t (Et+m t 
the decomposition of the squared scattering amplitude |.M| 2 can be written as 



|M| 2 = a + ^ S ? + &^ + c^ St ^, (20) 

and by comparing expressions ( [T8| ) and (201 one can extract the functions A, £?* {B\) and Cy. With this at hand, 
the various top spin observables (O) can be calculated as 

{0(S t , S f ))z = / dr tf Tr[ P/ • 0(S t , Sf)] , (21) 

where 07 = J" dT ( jTr[/?f] is the unpolarized production cross-section and St = er/2 (8) 1 (S^ = 1 ® c/2) is the top 
(anti top) spin operator. In particular, we consider the following spin observables 

01 = St-Sj, 

2 = S t a, 2 = S f -b, 

3 = 4(S t -a)(S r b), (22) 

which give the net spin polarization of the top - anti top system, polarization of the (anti) top quark, and the top - 
anti top spin correlation both with respect to spin quantization axes a and b, respectively. At the parton level O3 is 
related to the spin correlation function CV, in Eq. (18 1, namely 

, n ] g«-(tt) + vtt(ll) - ggCNj - gjgCjjO 

[ 3> train) + <rti(U) + <ru(U) + vattft ' 1 J 

where the arrows refer to the up and down spin orientations of the top and the anti top quark with respect to the 
a and b quantization axes, respectively. It can be measured using the double differential angular distribution of the 
top and anti top quark decay products: 

-j 7TT 7T = 7 (l + 5 t cos0/+Bjroos0 f -Ccos0/cos0f) , (24) 

a dcosOfdcosOf 4 v 1 11 

where 9f(0f) is the angle between the direction of the top (anti top) spin analyzer /, (/) (which can be either a direct 
t (t) daughter W + ,b (W~ ,b) or a W /+ (W /_ ) decay product i + {i~),v{v) or jets) in the t (t) rest frame and the a 
(b) direction in the ti center of mass frame (when the corresponding frame transformation is a rotation free boost, 
c.f. [H]). Analogously O2 and 62 are related to the (anti)top spin polarization coefficients B t and B^. We note in 
passing that in absence of CP violation B = B t = TBi for a = ±b. For perfect (anti)top spin analyzers whose flight 
directions are 100% correlated with the directions of the (anti)top spin then 

(0 3 >=C, (0») = B t , {0 2 )=B- t . (25) 

This limit is a good approximation for the charged leptons from W decays 53J. For other (anti)top spin analyzers 
one needs to apply the corresponding top spin analyzing power factors Kfif\ (where Ki+(i~) — 1) in Eq. (24 1 as 

C^CKfKf, B t (i) -> B t(i)K f{ fy (26) 

The values of are presently known at NLO in QCD and can be found in |53j . Finally, 0\ can be probed using 
the spin analyzer opening angle distribution 

- ~T~1 = ^ (! ~ ^ cos ^) 1 (27) 
a a cos 2 
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where <j> is the angle between the direction of flight of the two (top and anti top) spin analyzers, defined in the t and 
t frames, respectively. Again, for / = £ + , f — £~ one then obtains 

(Oi)=D, (28) 

while for other spin analyzers the appropriate Kf/f\ corrections should be applied. 

The arbitrary unit vectors a. and b specify different spin quantization axes which can be chosen to maximize the 
desired polarization and correlation effects. We work with the following choices: 

a = b = ki , ("helicity" basis) , 

a = b = p , ( "beamline" basis) , 

a = b = dx , ( "off — diagonal" basis, specific for model X) , (29) 

where p is the direction of the incoming beam and ki is the direction of the outgoing top quark, both in the tt center 
of mass frame. 

By the detailed study of the top (anti top) decay products one can obtain valuable information about top spin 
observables and use them to distinguish among the different models addressing the FBA puzzle. In order to maximize 
top spin effects it is advisable to choose a proper spin quantization axis. For the leading order QCD quark-antiquark 
annihilation dominating the tt production at the Tevatron, a special off-diagonal axis was shown to exist |54| . for 
which the top spins are 100% correlated. It is given by quantizing the spins with the axis dg^sM determined as 

h _ ~P + (l-7)zki 

d ^ SM - Vl-(l-7 2 )^ ' ( } 



where z — p • ki = cos 6 and 7 = E t /m t = l/^/l — j3 2 and interpolates between the beamline basis at the threshold 
(7 — > 1) and the helicity basis for ultrarelativistic energies (7 — > 00). There is no such optimal axis for the gluon-gluon 
fusion process (dominating the tt production at the LHC) 5 , but it is always possible to find a basis, in which spin 
correlations are maximal. A general procedure for finding such an off-diagonal basis is given in [561 15 7) . The idea is 



to determine the maximal eigenvalue of the matrix function Cjj in Eq. ( 18 ) and the corresponding eigenvector, which 
provides the off-diagonal quantization axis dx- Such an off-diagonal axis can be constructed for any model. As an 
example we present the explicit formula for the axigluon model (see also [5]) 

, -P+((l-7) Z -d^^) ki 

O-A = j , (31) 

y/l-(l- l 2 )z 2 + 2 1 ^^Tg A s /(s - m\)z - (1 - 7 2 )</V7(s - O 2 

where 9 a and vtlq are free parameters of the axis. The qq annihilation induced top spin correlations in this basis 
are maximal when these parameters match the actual axigluon resonance parameters (i.e. when qa = 9a9a an d 
771^4 = tog). It turns out however, that such custom axes are not necessarily better in discriminating between the 
SM and NP contributions to top pair production via qq annihilation. This can be easily understood by treating the 
SM and SM+NP as two competing hypotheses, among which we wish do discriminate using top spin correlations. 
Then the SM qq off-diagonal axis is already optimal for one of the hypotheses (it maximizes the leading QCD qq 
induced spin correlations) and is furthermore completely fixed by known SM parameters. Conversely any NP model 
off-diagonal axis needs to reduce to the SM one in the limit where the NP parameters determining the NP off-diagonal 



basis decouple (gA — ► and/or to a — > 00 for the axigluon axis in Eq. (31)). In the case of the axigluon model off- 
diagonal axis, we have checked explicitly, that in the interesting regions of parameter space and for arbitrary values of 
9a, fnA, the axigluon predictions do not deviate from SM values significantly more than in the helicity or the SM qq 
off-diagonal axis. Consequently, in our numerical study we use the helicity (for spin correlation Choi and polarization 
Bhei coefficients), beamline (for Cbeam, -Bbcam) and SM qq off-diagonal (for C ff, B q) spin quantization axes defined 



in Eqs. (29) and (30) respectively. 



At low m t i the top quark pair production via gluon-gluon fusion is dominated by like-helicity gluons. Consequently, spin correlations 
are maximal in the helicity basis |55) . 
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FIG. 7: Correlations between the NP contributions to the inclusive FBA and various spin observables at the Tevatron (see text 
for details and definitions). The present experimental results (68% C.L. regions) are shaded in horizontal and vertical bands. 
The NP model predictions are determined from the global fit as specified in Sec. IV and are bounded by full (axigluon G' in 
the low (mo < 450 GeV in black) and high (jug > 700 GeV in gray) mass regions), dashed (scalar color triplet A), dotted 
(scalar color sextet E) and dot-dashed (neutral component of the scalar isodoublet cj>° in the low (m^ < mt in darker shade) 
and high (m^ > 200 GeV in lighter shade) mass region) contours. 



B. Results 



In this section we present predictions for the various top spin observables at the Tevatron as well as the 7 TeV (and 
8 TeV) LHC within the various NP model parameter regions which are able to address the FBA puzzle, as determined 
in Sec. |IV| In particular we present correlations between the inclusive and high m t i FBA values as measured at the 
Tevatron, and the shifts of the various spin observables from their corresponding SM values. We define (see Sec. Ill) 
AA FB = A FB - Afg, AC t = Ci- Cf M and AD = D - D SM . On the other hand since QCD produced top quarks are 
not polarized, (neglecting tiny electroweak contributions) we assume Bf M ~ and present results for Bi in presence 
of NP directly. The predictions for the relevant spin observables at the Tevatron are shown in Fig. [7] First note 
that the results for the SM qq off-diagonal axis at the Tevatron turn out to be almost identical to the beamline axis 
(and very similar at the LHC, see Fig. [8]). Both bases provide good potential discrimination between color sextet on 
one hand, and color triplet or isodoublet scalar models on the other hand. The off-diagonal basis exhibits marginally 
better sensitivity only for the axigluon (G") model. However, since purely axial couplings of G to quarks do not 
produce polarized top quarks, Bi vanishes for the axigluon model and consequently we do not plot B a fi dependence 
separately. 

We observe that existing spin observable measurements at the Tevatron do not overly constrain selected NP mod- 
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FIG. 8: Correlations between the NP contributions to the inclusive FBA at the Tevatron and various spin observables at the 
7 TeV LHC (see text for details and definitions). The present experimental results (68% C.L. regions) are shaded in horizontal 
and vertical bands. For AChei we also show the 95% C.L. contour in thin dashed line. The NP model predictions are determined 
from the global fit as specified in Sec. |IV| and are bounded by full (axigluon G' in the low (mc < 450 GeV in black) and high 
{tug > 700 GeV in gray) mass regions), dashed (scalar color triplet A), dotted (scalar color sextet E) and dot-dashed (neutral 
component of the scalar isodoublet <j>° in the low (m<j, < m t in darker shade) and high (m^, > 200 GeV in lighter shade) mass 
region) contours. 



els. Some sensitivity to the light scalar isodoublet model is exhibited by the recent beamline axis spin correlation 
measurement by D0 as seen in the center left plot in Fig. [7] On the other hand (anti)top polarization 
both in the beamline and in the helicity basis) offers a very powerful probe of scalar i-channel models and a O(20%) 
precision measurement (in helicity basis) could already test (and discriminate between) the scalar color triplet (A) 
and isodoublet (<^°) model explanations of the FBA. Finally, the axigluon (£?') models in general give very small 
contributions to the chosen spin observables. For example, at the Tevatron, spin correlation measurements at 0(2%) 
precision would be required to probe such FBA explanations. 

The results for the relevant spin observables at the 7 TeV LHC are shown in Fig. [8j 6 Among these, presently the 
most powerful probe of FBA inspired models is the helicity basis spin correlation as measured recently by ATLAS [U] . 
In particular it already represents a non-trivial constraint for the scalar isodoublet and heavy axigluon models. In 
the light scalar isodoublet scenario, the large negative deviation in AChoi can be traced to sizable non-standard 



6 The results for AD, Ad and Bi at the 7 TeV and 8 TeV LHC are almost identical and we do not show the later separately. 
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t — > cj> u decay rate, reducing the experimental tt reconstruction efficiency (which is then compensated by enhancing 
the uu — > tt cross-section contribution). In term this leads effectively to a reduced gluon fusion component in the tt 
cross-section which in this basis contributes sizable positive spin correlations in the SM. 

Comparably sensitive (and presently unmeasured) observables are also the opening angle spin correlation coefficient 
D, beamline axis spin correlations (Cb oa m) and hclicity axis top polarization (_Bhoi)- 0(5%) precision measurements 
of these quantities could discern among the i-channel scalar models. Finally, again we observe that the axigluon (C) 
model (especially in the light nig < 450 GeV region) gives very small contributions to the chosen spin observables at 
the LHC and will be difficult to probe in this way. 

VI. CONCLUSIONS 

We have performed a comprehensive analysis of tt production phenomenology at the Tevatron and the LHC within 
effective single NP resonance models addressing the FBA puzzle. We have quantified an inherent tension between the 
large positive FBA measurements at the Tevatron and precise CA measurements at the LHC (consistent with zero). 
In particular the later conclusively exclude the W' and Z' explanations of the FBA anomaly. Among the considered 
models, only a color octet axial-vector axigluon state (of mass m,Q ~ 400 GeV or tug ^ 1 TeV) can reproduce the 
central experimentally determined values of the inclusive and the high m t t bin FBA without being in severe conflict 
with other constraints. In addition, a light scalar isodoublct (m^ < m t ) model predictions can barely reach the one 
sigma region for the inclusive FBA average, while the central value is in conflict with the LHC cross-section and CA 
measurements and also the recent top spin correlation measurement by ATLAS. Finally, the scalar color triplet and 
sextet models are constrained by cross-section and m t i spectrum measurements at the Tevatron. However, given the 
caveats associated with properly evaluating the contributions of these models (and also the isodoublet model) to the 
existing experimentally reconstructed m t i spectra (unfolded using SM signal templates), we suggest the experiments 
at the Tevatron and the LHC perform dedicated studies to settle the issue. 

For the favored NP model parameter regions we have derived predictions of (anti)top spin polarization and tt spin 
correlation observables at the Tevatron and the LHC. All scalar models addressing the FBA puzzle predict significant 
deviations in several top spin observables. At the Tevatron the most promising are the (anti)top polarization fractions 
with respect to the helicity spin quantization axis, which can deviate by more than 20% from their SM predicted 
values and their measurement could discriminate between scalar isodoublet and color triplet models. At the 7 (and 
8) TeV LHC the scalar t — (u— )channel models predict smaller effects, but helicity axis top spin correlation and 
polarization measurements at the 5 — 10% accuracy can yield competitive constraints, as exemplified by the recent 
ATLAS [41] result. On the other hand, a light [juq < 450 GeV) axigluon model predicts very small effects in top spin 
observables both at the Tevatron and at the LHC, and will be difficult to constrain in this way (much larger effects 
may be expected only if in addition to the axial component, a vector coupling to qq pairs is present [51 151j). 

While in the present top spin observables' analysis we have focused on inclusive tt samples, binned distributions of 
these quantities in ra t i and (anti)top rapidities could further enhance NP effects relative to the SM dominant gg — > tt 
subprocess and help to discriminate between various possible contributions (c.f. [5J [SHU 1581 159] for recent proposals 
along these lines). 
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Appendix A: Partonic Amplitudes 
1. Parametrization of Standard Model and New Physics Contributions 

At the partonic level in the SM the q{p\)q{p2) — >• £(fci, St)^(^2, Sf) amplitude is given by 

iMf M (s t ,St) = i^v( P 2WT a ufa)u(k 1 ,at)i li T a v(h 2 ,8 { ) , (Al) 

with s = (pi +P2) 2 — (fci + k^) 2 , t = (px — ki) 2 and u = (pi — fe) 2 . Also, g s denotes the QCD coupling while T a are 
the Gell-Mann SU(3) generator matrices. 
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Then the unpolarized qq SM cross section is 



da 



gg,unpol. 

SM 

dt 



2na 2 s 



[2- (1 - cos 2 9)/3 2 ] 



while for gg initial state 

ua SM 

dt 



48s 2 (1 - P 2 cos 2 9) 2 



p 4 cos 4 9} 



(A2) 



(A3) 



where, as usual, f3 = y 1 — 4m| / s and 9 is the scattering angle in the ti center of mass system. 

Since we are also interested in top spin observables, we present also the results for the polarized cross sections. The 
new physics models which we consider have particles which couple exclusively to quarks and antiquarks and therefore 
their amplitudes interfere only with the SM qq parts. The cross sections receive the contribution from the SM, the 
interference between the SM and the NP, and the NP part only. For all models, denoted by X, we thus write 



J 99,Pol. 
aa SM+X 



1 



\M 



pol. 



( pol.|2 _ a 



where N — 3 is the number of colors each square of the polarized amplitudes can be written as 
and further 

b^l = b\{ Pl - S t )+b 2 (p 2 - s t ) 

&M*f = b{{pi-st) + bl{p2-st) 



c(s t ■ st) + d[(jpi ■ s t ){pi ■ + (pa • s t )(p2 ■ s^} 
+e(pi • s t )(p 2 ■ st) + f(p 2 • s t )(pi ■ St) . 



Due to the CP invariance in our models, b\ 



-b\ and b\ 



The spin coefficients for the polarized qq squared amplitude in the SM are as follows: 

<i<l 



«SM = Cf M [2-P 2 (l-z 2 )] 

4V = cf M [p 2 {i-z 2 )] , 

4 

-( 
s 

4 



J SM 



f~,qq 
^SM 



f<?<? _ 
'SM — 



C 



'I 'I 
SM 



■~(1 + Pz) 
-(1-/3*) 



(A4) 

(A5) 

(A6) 
(A7) 

(A8) 



(A9) 



where Cg 9 M — 167r 2 a 2 (7V 2 — 1) and z = cos 9. All other coefficients are zero. On the other hand, the polarized SM 
<7<7-cross section has the form 



da 



gg, pol. 

SM 



dt 



16tts 2 4(iV 2 - 1) 



and |Mg M pol '| 2 can be written with the help of 



a gg _ 
-sm — 



J SM 



r<gg 



r>gg 



-/? 2 (l + /?z)(l-z 2 ) 



/3 2 (l-/?z)(l-z 2 ) 



(A10) 



4m = Cg 9 M [l + 2/3 2 (l~z 2 )-(2-2z 2 +z 4 )/3 4 
4 9 M = Cf M [l-2/3 2 + (2-2z 2 + z 4 )/3 4 ] , 
4 



(All) 



where C s q M = 327r 2 a 2 ^ N //jT^^p ^ ^ N N ^ ■ These expressions can be also found in Sees. 2.1.1. and 2.1.2 of 
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2. New Physics Models 

a. Axigluon 



The relevant interaction part of the Lagrangian describing new physics in the tt production due to s-channel axigluon 
G' exchange is given by , 



= -Q(9 9 v ~ 9 g A l5)€'q - t(g v - g Al5 )G't , 



The new qq — > tt amplitude is then 



s — rriQ + imc^G 



v\p2)l^{g v +g\l5)u 1 {p 1 )u\k ll s t ) 1 ^g v +g\ 1 r )v\k 2 ,st) 



(A12) 



(A13) 



where mg is the axigluon mass and Tq is the manually introduced axigluon width. The resulting unpolarized qq — > tt 
cross section is of the form 160 



da 



gq,unpol. 

SM+G' 

dt 



dt ' + 9s 2 a ° I (s - m 2 ,) 2 + T 2 G m 2 



S( " '" r>; ' 2 (g y g v (l + /3 2 z 2 + 4m 2 / S ) + 2g\g\Pz) 



G G 

(g _ m2 ;f +r2m2 [((^) 2 + (3%?) ((Str) a (l + £V + 4m?/ S ) 

V G / G G 

+(^) 2 (1 + /3 2 z 2 - 4m 2 / S ) + 8/3*^^)] 1 , 



(A14) 



where now g$\ = g'yjg. 



The polarized parts are given (for g v = g\r = 0) by the following expressions, eq. (A8) 

a SM*G' — Csm*G' [4/3 , 
eSM*G' = Csm*G 



and 



where 



16m 2 



SM*G' 



— Csm*C 



-16m 2 



a G > 

C G ' 

e-c = C G 

fc> = Cg 



Other coefficients in ( A8 1 are vanishing. 



C G '[l-^m 2 / S + f3 2 z 2 }, 
C G 4-l + 4m 2 /s + (3 2 z 2 
4 
s 
4 
s 



(1 - 4m 2 /s + I3z) 
(-1 + 4m 2 /s + f3z) 



r< ia^ 2 „ 2 (AT 2 i n 9a9a s ( s ~ m G) 



C G > = 16ttX(A^ - i) 



(&) 2 (&) 2 « 2 

(s — m? G ) 2 + fgnig 



(A15) 



(A16) 



(A17) 



b. Z' and W 



The interactions of the heavy weak bosons Z' and W contributing in the t-channel of the tt production are given 
by the following Lagrangians: 



/»int. 
>-7' 



-ujME p L + fiP R )tz' + h.c, 



(A18) 
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The relevant amplitude induced by the new physics contributions is 



(A19) 



iM z > = iu(k u s t )^(f£ 'P L + fg P R )u( Pl ) [ - 



Z' ; 



t — m z , mr z , (t — m z , ) 



v(p 2 ) 7 v (fE'p L + fi'p R )v(k 2 , Si ) 



(A20) 



where t M = (pi — ki)^ , and similarly for the interaction with the W, by the exchange Z' — > W and u — » d. Then 



j uu,unpol. 
m SM+Z' 

dt 



da 



im.unpol. 

SM 

dt 



+ <*s_ (fE ) 2 + (/D 2 {2m 2 + 2gm 2 + + m 2 s)} 



9s 3 i — m\i 



+ 



16tts 2 (t - 



^ 2 { ((/f ) 4 + (/I') 4 ) «? + 2(/f ) 2 (/I?.s( S - 2m?) 



+ 



rag 
4m| 



;((/f) 2 + (/I') 2 ) (t 2 + 4 S m : 



(A21) 



where u t = s(l + /3z)/2,t t = s(l - /3z)/2. 



The polarized case we consider only the right-handed couplings and the non- vanishing spin coefficients in ( A8 ) are 

™2 



a-SM*z' — Csm*z> 
csm*z' — CsM*Z' 
zsm*z' — Csm*z' 



4to 2 + «(1 + /^) 2 + ^H 4 ™? + s (! - Z 32 ) 2 ) 



-4m? + s(l - /3 2 z 2 ) - ^J;(4m? - ,(1 - /3 2 z 2 )) 
- f2m 2 - S (l + /3z)-^-^(2m 2 + S (l + /3z)) 



SM*Z' 



C. 



SM*Z' 



2m 2 t + s(l - (3z) - ^(2m 2 - s(l - f3z)) 



(b\)sM*Z' — Csm*Z' 



SM*Z 



> = C, 



SM*Z' 



and 



az> = C z > 
cz> = C z > 



-■2m, ( l + /3z+ ^_|-(3-/3z) 



-2m t (^3 + /3z+^-^(l-/3z) 



,s(l + / 3z) 2 + 4-^ + -^ S (l-/?z) 



^( S (l-/3V)-4m t 2 ) 



(bi)z' 



Cz> 
C z > 
C z > 



m z , 
4 m? 



2m 2 - s(l + 0z) 



TO 2 , 



l^( 2m 2_ s(1 _^)_ 4ro 2 ) 
S TO^/ 



-2m t - 



mr 



mr 7 . 



11)7 



(l + pz) + ^(l-Pz) 



11)7 



--!'»/ ( H - ,:, + 2%( 1 -' 8z ) 



with 



Csm*z' = 87ra s 



£ - m|, 2 



(m|, - ty 



N 



The W' case is again given by the same expressions by substituting Z' — > everywhere. 



(A22) 



(A23) 
(A24) 
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c. Scalar isodoublet 

The interactions of the weak doublet scalar $ ~ (1,2) 1/2 with quarks are given by [TBI [61] 

= -VijQLiURj^ - /;;>// ,</;,■, <i- + h.c. , (A25) 

where $ = 2x2$*. The amplitude coming from the exchange of the neutral isodoublet component (jP is then (with 

Vij = Vij) 

lM * = -^ J^\^ a {k u s t )^{l~ l5 )v\k 2l s l )v\p 2 ) 1 ^l + l5 )u a {p 1 ) 1 (A26) 
where m$ is the 4>° mass, resulting in the cross section 

1 tiSjUiipol. tifi.unpol. I 12 2 i/2 I 14 i / 2 j-\2 

d(J SM+4> _ da SM p a s \y ls Y m z t s + {m z t -t) z I2/13 1 1 - *) 



di df 9 s 3 - < 64tt s 2 (m 2 - 1) 5 

The coefficients in the polarized cross section are 



and 



where 



1 6 TOs |y 13 | 2 N 2 -l |2/i 3 | 4 



(A27) 



asM*<p = C SM *d> [t 2 + m 2 (s - 2t) + mf] , 

csm*4> = Csm*4> [ lu — mt] j 

esM*<f> — Csm*4> [(— 2)(s + 1)] , 

fsM*4> = Csm*4> [2t] , 

(&*i)sm*0 = Csm*4> [m t (u + 2t - 3m 2 )] , 

(6 2 )sm*0 = C 5A f*0 [mt(t- m 2 )] , (A28) 



a<?i = [(m 2 - t) 2 ] , 
e = C*0 [-4to 2 ] , 
[b% = C [2m t (t-m 2 )] , (A29) 



Csm*4> — — ( + !Jli ^2\ — o — ' ^4> ~ 7Z^T^7y2 ^ 2 ' (A30) 
The rest of the coefficients are vanishing. 



d. Scalar color triplet and sextet 

The relevant interaction Lagrangians in this case are given by 

4f = -g iA)ij e abc u a Rti (u b RJ ) c A c + h.c. , (A31) 
with A c being a (3, 1, —4/3) state, and 

4 nt - = (4,(4,,)° + ^(^) C )S oi>t + h.c , (A32) 

with Yj ah being a (6, 1,4/3) state. For the color triplet amplitude one finds 

t M A = i 1 ^^ -^[u a (k 1:St )r(l + l 5 )v a (k 2 ,s i )v b (p 2 ) 1 ^l + l5 )u b (p 1 ) 

-z2 a (fc 1; St )7 M (l + 75)^(^2, st)v»(p 2 ) lfl (l + j 5 )u a ( Pl )} , (A33) 
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while for the color sextet case we have 



(A36) 



. Ms = . lg(E)13l 1 fcafo ) ^(l + TaXte^fo) 7M (l+7 5 y( Pl ) 

8 u — 

+w a (fei, s t )7 M (l + I5)v b (k 2 ,st)v b (p 2 ) lfl (l + 7sK(pi)] • (A34) 
The cross section for the color triplet scalar is then given by 

^ unpoL a fl |g(A)i 3 | 2 m 2 s + (mf - u) 2 |g (A)13 | 4 1 (mf - u) 2 
dt dt 9 s 3 to^ - u 48tt s 2 ( m 2 _ u )2 " ^ J 

In the case of the color sextet model one finds instead 

^S+T L = dtT^r P °'- | a s \9^u\ 2 m 2 s+(m 2 -u) 2 | |g (s)13 | 4 1 (to 2 - u) 2 
dt dt 9 s 3 m| - it 24-7T s 2 (m| - u) 2 ' 

The nonvanishing polarization coefficients for these models are: 

Osa/*a/s = C S m*A/T. [u 2 + m 2 {s~2u)+mf\ , 

c SM*A/£ = C,SM*A/£ [t u ~ m t] ) 

esM* a/e = Csm*a/s [2tt] , 

/SM*A/E = C S M*A/S [— 2(s + li)] , 

(61) sM*A/E = C , SM *A/E ["l f (u - TO 2 )] , 

(6 2 ) sa/*a/s = C SMtA/s [±m t (3m 2 -t-2 W )] , (A37) 

and 

«A/s = C A/S [(to 2 - u) 2 ] , 
eA/s = C A/S [-4to 2 ] , 
(6*) A/S = ±C A/S [2to 4 (to 2 - u )] , (A38) 

where 

Csm*a/ S - ± - (u _ m2) (JV - 1) , C A/E - u)2 [2N(N T l)\ ■ (A39) 
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